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Low light level charge coupled devices (L3CCDs) have recently been developed, 
incorporating on-chip gain. They may be operated to give an effective readout noise 
much less than one electron by implementing an on-chip gain process allowing the 
detection of individual photons. However, the gain mechanism is stochastic and so 
introduces significant extra noise into the system. In this paper we examine how best 
to process the output signal from an L3CCD so as to minimize the contribution of 
stochastic noise, while still maintaining photometric accuracy. 

We achieve this by optimising a transfer function which translates the digitised 
output signal levels from the L3CCD into a value approximating the photon input as 
closely as possible by applying thresholding techniques. We identify several threshold- 
ing strategies and quantify their impact on photon counting accuracy and effective 
signal-to-noise. 

We find that it is possible to eliminate the noise introduced by the gain process at 
the lowest light levels. Reduced improvements are achieved as the light level increases 
up to about twenty photons per pixel and above this there is negligible improvement. 
Operating LSCCDs at very high speeds will keep the photon flux low, giving the best 
improvements in signal-to-noise ratio. 

Key words: instrumentation: detectors - techniques: photometric - methods: sta- 
tistical - methods: numerical. 
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1 INTRODUCTION 

Charge coupled devices (CCDs) are ideal detectors in 
many astronomical applications. They are available in large- 
format arrays, have high quantum efficiency (QE), a linear 
response, and, if cooled sufficiently, a very low darf; cur- 
rent. Their major shortcoming is readout noise, i.e. the addi- 
tional noise added by the on-chip output amplifier, where the 
charge of the detected photo-electrons is converted into an 
output voltage. Currently, the typical noise levels achieved 
at slow readout rates (e.g. kHz pixel rates) are li ttle bet- 
ter than ~ 2 electrons per read jjerram et al ]|200l'). At the 
higher readout rates (MHz pixel rates) often used, for ex- 
ample, for adaptive optics and interferometric applications, 
far poorer performance is the norm, with typical noise le vels 
of ~ 10 — 100 electrons per readout dJerram et al.ll200ll) . 

* E-mail: abasden@mrao.cam.ac.uk 



A novel solution to this problem, in which on-chip gain 
is used to amplify the signal prior to readout, has recently 
been demonstrated by E2V Te chnolog i es (fo rmerly Mar- 
coni Applied Technologies. I Jerram et al. i200ll) ). In this ap- 
proach, an extended serial register is used to allow electron 
avalanche multiplication so that a large mean gain can be 
realised prior to a conventional readout amplifier. Although 
the effective gain can be very large, the detailed process by 
which the signal is amplified is stochastic and so introduces 
additional noise at the output. The effects of this noise, and 
its correction by a judicious analysis of the output from a 
typical low light level CCD (hereafter L3CCD) are the pri- 
mary subjects of this paper. 



The content of the paper is as follows. In Section 2 we 
describe how the L3CCD works and develop a numerical and 
probabilistic model of the gain mechanism. In Section 3 we 
discuss techniques for analysing the noisy output signal of 
an L3CCD and the implications for photometric accuracy 
and signal-to- noise ratio (SNR) , and we conclude in Section 
4. 
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Figure 1. A schematic diagram of an L3CCD showing the mul- 
tiphcation register. 



2 CHARACTERISTICS OF THE L3CCD 

2.1 Principles of L3CCD Operation 

The L3CCD architecture is similar to that of a normal CCD 
except that it has extended serial register, called the mul- 
tiplication register, allowing for additional serial transfers 
before the signal reaches the ouput amplifier (see Fig. 0. 
The electrode voltages in the multiplication register can be 
adjusted so that avalanche multiplication of the electrons 
occurs as they are moved through each element of this part 
of the device. At each step, the probability, p, of producing 
an additional electron per initial input electron is small — 
typically this may be 0.01 — 0.02 — but the cumulative effect 
of many transfers can be very large. For example, for a reg- 
ister comprising r = 591 elements (the number of elements 
in the CCD65 from E2V Technologies), the mean gain, g, 
will be (1-l-p)'' = 6629, when p = 0.015. 

In general, the probability distribution of the output x, 
for an input of n (integer) photo-electrons, where the mean 
gain is g will be given approximately by (Appendix A and 
Fig.HJ 



' ^exp(-x/g) 



(1) 



when the photon input level is relatively small and the gain 
is large. This distribution has a mean of ng and a variance 
of ng^ and at high light levels is approximately Gaussian. 
The signal noise introduced by this multiplication process 
is independent of the photon input. The SNR of an L3CCD 
output is obtained by combining in quadrature the noise due 
to the Poisson nature of light with the noise from the multi- 
plication process. At a mean light level /i photons per pixel, 
this gives an SNR equal to y/Ji/ E where E (the excess noise 
factor, hereafter ENF) is equal to \/2 for an L3CCD output 
with large gain. The nature of the multiplication process 
means that in general there will not be a one-to-one map- 
ping between the number of electrons entering and leaving 
the multiplication register (see Fig. so that in principle 
there will be always be some uncertainty when estimating 
the input flux. For a CCD with no gain, E is equal to unity, 
and so to achieve the same SNR when using a large gain, we 
will need to detect twice as many photons, meaning that the 
multiplication process has effectively halved the QE of the 



0.00016 
0.00014 
0.00012 
0.0001 



6629 13258 19887 26516 33145 39774 46403 53032 59661 
Gain 

Figure 2. L3CCD output probability distribution for a given 
photon input. Inputs of 1 — 4 photons are shown for a mean gain of 
6629. Note the significant overlap between different inputs which 
causes uncertainties when estimating the input from the output. 



device. This is a serious loss as astronomers are happy to 
spend a lot of money on optical coatings to increase system 
throughput by only a few percent, and do not want to accept 
the loss in effective QE that the use of an L3CCD implies. To 
be able to increase the effective QE using signal processing 
is therefore essential if L3CCDs are to be used effectively 
for astronomy. This is the main driving force behind our 
investigation, which aims to use our additional knowledge 
of the system to allow us to reduce the uncertainty in the 
measured fluxes. 



2.2 Reducing the noise 

We consider processing strategies which allow up to mini- 
mize the stochastic multiplication noise. At very low light 
levels, much less than one detected photon per imaging pixel, 
we are able to use the L3CCD in a photon counting mode. 
At these low light levels, we will either get zero or one pho- 
tons in a pixel per integration time. This will then result in 
either a signal much smaller than the amplifier read noise 
(zero photons), or when a detected photo-electron is ampli- 
fied, a much larger signal. Provided the mean multiplication 
gain is much greater than the read noise, we can then treat 
any signal above some threshold as having arisen from one 
photon event. Replacing the range of signals we get from a 
pixel containing one photon (curve 1 in Fig. \^ eliminates 
the variance in output signal introduced by the multiplica- 
tion process of the L3CCD and increases the SNR to a value 
we expect from a conventional CCD with negligible read-out 
noise. 

At higher light levels this technique will not work, since 
coincidence losses (more than one photon falling on a pixel, 
being interpreted as one photon) will become increasingly 
important. Nevertheless, the idea of thresholding can still 
be helpful if we use more than one threshold. We can see 
how this might work as follows. An L3CCD with a mean 
multiplication gain of g might give an output signal of fg. 
This could be due to a single detected photon that has been 
amplified by an unusually large amount. Alternatively, if we 
set thresholds Tn and r„+i and estimate the flux to be n de- 
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tected photo-electrons before multiplication if the detected 
signal satisfies T„ < fg < Tn+i, then we will be able to 
maximize our chance of estimating the flux correctly if we 
choose thresholds correctly, ft is the determination of these 
threshold levels and any corrections that need to be applied 
in order to conserve flux which we now consider, together 
with their impact on the L3CCD output SNR. 



2.3 Thresholding schemes 

We have considered many different signal-processing strate- 
gies, and provide details of the most useful selection here. 
We investigate which of these provide the largest SNR im- 
provement. When investigating these strategies, we assume 
that the data is first thresholded at some level above the on- 
chip readout noise level (typically 6(t where g is the RMS 
noise due to the readout amplifier) so that amplifier read 
noise is negligible. Following this, the strategies considered 
are: 

1. Analogue: The output signal is divided by the mean 
gain. 

2. Photon Counting (PC): If the output signal is above a 
single fixed noise threshold it is treated as representing one 
input photon. 

3. Poisson Probability (PP): Threshold levels are set 
where the probability of an output resulting from a mean 
input of n Poisson photons is equal to the probability of 
the output resulting from a mean input of n + f Poisson 
photons. 

The Photon Counting (PC) thresholding strategy will 
underestimate flux at light levels where there is a non-zero 
probability of more than one photon being detected on a sin- 
gle pixel. The Poisson Probability (PP) strategy will over- 
estimate the flux, particularly at low light levels, and so 
a theoretically determined correction is necessary after the 
detections, which we investigate in section 3. 



2.3.1 Threshold boundaries 

At a given mean light level, ^ photons per pixel, the L3CCD 
output with mean gain g can be estimated by providing the 
photon input in Eq. with a Poisson probability distribu- 
tion, giving: 



Table 1. Threshold boundaries in units of mean gain, for a PP 
thresholding strategy. 



p(^^^)^^^^p(-M-x-/ff)M"(^/gr 



g{n — f )!n! 



(2) 



where p{x,fi) is the probability that the L3CCD output will 
be X when the mean light level is ^, and the mean gain is g. 

We use this distribution to determine threshold bound- 
aries for the PP thresholding strategy, which are given in 
Table 1. These are placed at the points where the probabil- 
ity of getting an output x with a mean light level /i is equal 
to the probability of getting the same output x with a mean 
light level of /^-l-f as shown in Fig.|H] i.e. finding x such that 



p{x,a) =p{x,a + l) 



(3) 



where the a' threshold boundary is placed at position x, 
and p{x,a) is defined in Eq. |5| This results in threshold 
boundary positions independent of the mean light level. The 
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Figure 3. The L3CCD output probability distribution for mean 
(Poisson) light levels of 1 (loft curve) to 4 (right curve) photons 
per pixel, as given by Eq. |^ Threshold boundaries for the PP 
thresholding strategy are placed where these curves cross. 



theoretically predicted mean flux and ENF arising from this 
strategy is also calculated. 



2.4 Threshold evaluation 

Monte-Carlo simulations were used to verify that our thresh- 
old positions were chosen correctly, and that our calculations 
of photon input estimation and ENF were correct. Input 
photon streams were generated for differing mean light lev- 
els assuming Poisson statistics, and were injected into the 
first pixel of a multiplication register of length r. The trans- 
fer of this signal to the next pixel of the register was then 
computed assuming a probability p that any given photo- 
electron would be amplified to give 2 photo-electrons, and a 
probability 1 ~ p that the transfer took place with no ampli- 
fication. This process was then repeated a further r — 1 times 
to simulate the output expected for the initial input signal. 
In this way, the expected output for L3CCDs with differing 
multiplication probabilities and register lengths could easily 
be generated for differing input light levels and threshold 
boundary positions. Large numbers of these output data se- 
quences were then processed to verify that the theoretical 
ENF and flux estimation calculations were correct. 



2.5 Figures of merit 

In order to assess the performance of different signal- 
processing strategies, the quality of the signal recovery was 
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quantified using both the ENF and the following misfit func- 
tion, M: 



M 



(4) 



where ni is the true input photon count and yi the value of 
the input signal estimated from the thresholded L3CCD out- 
put, Xi. Data sequences comprising many tens of thousands 
of signal values were generated so as to allow different meth- 
ods for generating the estimates {yi} from the raw L3CCD 
outputs {a:;i}. The minimization of this misfit function cor- 
responds to the best input prediction. 

An ideal detector will give a misfit of M — 0, and a 
variance equal to the variance of the input (Poisson) data, 
H = n, with the SNR equal to y/Jl. A non-ideal detector will 
have greater dispersion in the output, resulting in a reduced 
SNR of ^/E where E is the ENF. 



3 DISCUSSION OF IMPROVEMENTS 

When investigating our L3CCD output thresholding strate- 
gies, it is helpful to consider three different light level 
regimes. The first of these is when the mean light level is 
low, much less than one photon per pixel per readout. Sec- 
ondly, we consider intermediate light levels, with between 
about 0.5 — 20 photons per pixel per readout, and finally 
above this, high light levels are considered. This separation 
allows us to apply the different processing strategies able to 
maximize the SNR at each light level. Our results are inde- 
pendent of multiplication register lengths typically found on 
LSCCDs (greater than 100 elements), though a shorter reg- 
ister length will generally give a slightly lower ENF (Eq.|5J. 

As shown in Fig. 01 and Fig. |S] the on-chip amplifier 
readout noise can have a significant impact on the flux de- 
termination if the mean gain is not much larger than the 
noise threshold level (below which any signal is ignored), 
particularly at low light levels. By setting the noise thresh- 
old level at least 4ct above the mean on-chip readout noise 
(assuming readout noise with an RMS of a), and using a 
mean gain at least ten times greater than the noise thresh- 
old level, we are able to minimize the effect of the on-chip 
readout noise. As light level or gain are increased, the effect 
of readout noise is reduced. 



3.1 Low light levels 

As discussed in section 12.21 if the mean light level is low 
(much less than one photon per pixel per readout), we can 
use a PC thresholding strategy, treating every signal above 
some noise threshold as representing one photon. This re- 
moves all dispersion introduced by the multiplication pro- 
cess and effectively eliminates any additional noise. The 
signal-to-noise ratio from many such samples scales as -y/n 
as shown in Fig.|^ The mean gain does not need to be accu- 
rately determined, but should be kept well above the readout 
noise, as mentioned in the previous section. If this is not the 
case then some real signals will have insufficient amplifica- 
tion and will be treated as noise, leading to inaccuracies in 
flux estimation. 

Coincidence losses will result in nonlinearities in the flux 
prediction as the light level increases. If the mean light level 




Figure 4. A graph showing the effect of on-chip readout noise for 
the PP thresholding strategy as a function of mean gain, for dif- 
ferent noise threshold levels. The upper curve has a noise thresh- 
old set at 1 X (T above the noise level, a = 50e~, a being the 
RMS noise, while lower curves set the noise threshold at 2a, 3(t, 
4(T (dotted curve), 5cr and 6a above the mean noise. A noiseless 
readout {a = 0) would give a value of unity for all gains. 
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Figure 5. A graph showing the effect of on-chip readout noise 
for a PP thresholding strategy when the noise threshold is set to 
4a above the mean noise level when a, the RMS noise, is equal to 
50e~. The upper (solid) curve is for 10 photons per pixel, lower 
curves for 1 (dashed), 0.1 (dotted) and 0.01 (dash-dotted) photons 
per pixel respectively. The effect of readout noise is minimized at 
higher light levels and gains. Curves are normalised to a value of 
unity representing a noiseless readout. 



is 0.2 photons per pixel then two or more photo-electrons 
will be detected on a pixel less than 2 per cent of the time, 
resulting in a relatively small coincidence loss which can be 
corrected easily. However, if the mean light level is 1 photon 
per pixel, coincidence losses are larger and we would esti- 
mate the light level to be only 0.63 photons per pixel. This 
nonlinearity can be determined and so we can correct the de- 
tected photon flux while this nonlinearity remains small, for 
light levels up to about 1 photon per pixel. At higher light 
levels, our estimated light level tends towards unity and so 
we are unable to deduce the correct light level without much 
uncertainty. 

Our PC thresholding strategy may be applied here, but 
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Figure 6. Excess noise factor as a function of light for the PP 
(soHd curve), PC (dotted) and analogue (dash-dotted) thresh- 
olding strategies. Curves are theoretical, while markers are from 
Monte-Carlo simulation, and negligible readout noise is assumed. 
An excess noise factor of \/2 is equivalent to halving the QE of 
the L3CCD. 
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Figure 7. Ratio of estimated to true input for L3CCD output 
data using the PP (solid curve), PC (dotted), corrected PP (using 
Eq. |5| dashed) and analogue (dash-dotted) thresholding strate- 
gies. Curves are theoretical, while markers are from Monte-Carlo 
simulation and negligible readout noise is assumed (g 2> O"). We 
see that flux is overestimated at low light levels using a PP strat- 
egy, but that this nonlinearity can be corrected for. 



will overestimate the flux, since there is a significant proba- 
bility that the output from a single photon input will be 
placed into the second (or greater) threshold bin. Fig. Q 
shows the size of the error in flux estimation for our PC 
and PP thresholding strategies as a function of light level, 
and hence the nonlinearity correction that should be applied 
after data from many frames has been thresholded in this 
way. 



3.2 Multiple thresholding at intermediate light 
levels 

At intermediate light levels up to about twenty photons per 
pixel, we cannot use a PC thresholding strategy as coinci- 



dence losses become large. However, as discussed previously 
it is still possible to process the output, reducing the ENF. 
Our PP processing strategy can be applied at any light level 
though is most advantageous in this light level regime, and 
gives decreasing improvements up to about 20 photons per 
pixel per frame. Above this, the photometric accuracy and 
ENF are indistinguishable from those obtained using the 
analogue processing strategy. 

Threshold boundaries from the PP processing strategy 
are independent of light level. Summing thresholded output 
signal values and applying a nonlinearity correction for light 
level provides us with a good estimate of the flux. Fig. Q 
shows the nonlinearity correction which should be applied. 
Without this, the flux will be overestimated at lower light 
levels since there is always a signiflcant probability that the 
output from a single photon will be interpreted as two or 
more photons. 

We flnd that theoretical results and those from our 
Monte-Carlo calculations agree almost perfectly, as would 
be expected, as shown for example by Fig. |7| 

3.3 High light levels 

At high light levels, the input photon (Poisson) distribution 
becomes symmetrical about the mean light level, having the 
form of a Gaussian. The multiplication noise distribution 
also tends to a Gaussian, and so the L3CCD output distri- 
bution is Gaussian. 

We treat the output as we would a conventional CCD, 
simply dividing by the mean gain, using the analogue pro- 
cessing strategy. This does not remove any of the dispersion 
int roduced by the mult iplication process, giving an ENF, E, 
of iMatsuo et al.lll985l) 

' ' (5) 
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l/r 



jH-l/r 



where r is the multiplication register length, g the mean 
gain, and the approximation is valid when r and g are large, 
as is usual for an L3CCD. The SNR is then yj'n/2, effectively 
halving the QE, and the misfit function (Eq. (|1J) is also 
greater than for other thresholding strategies, as shown in 
Fig.H 

This analogue processing technique may be applied at 
any light level, though due to the large ENF, other methods 
can give an improvement in SNR, particularly at light lev- 
els below about 20 photons per pixel per readout. At very 
high light levels, where the mean light level is a few times 
the square of the on-chip amplifier re ad noise of t he L3CCD 
(typically 50— 100e~ for fast readout . I Jerram ~ al. (2001)), 
the multiplication gain can be turned off, and then no sta- 
tistical noise is added, giving Poisson shot noise scaling as 
v'n. This is the mode in which a conventional CCD is used. 

3.4 Excess noise factors 

An ENF of ^2 effectively halves the QE of the CCD, which 
is a serious loss. Fig.|S|compares the ENF for various thresh- 
olding modes, as a function of light level. We can see that 
a combination of PC and PP thresholding strategies allows 
us to reduce the ENF for light levels up to about twenty 
photons per pixel. 
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Figure 8. Misfit as a function of Ugiit for the PP (solid curve), 
PC (dotted) and analogue (dash-dotted) thresholding strategies. 
Negligible readout noise is assumed. We see that a PC strategy 
should not be used at light levels greater than about one photon 
per pixel. pics2/misfitvlight.eps 



It is interesting to note that the ENF tends towards 
zero for the PC strategy at high hght levels. This is because 
we always interpret every signal as one photon, and so the 
result has no noise. However, in this mode, we are unable 
to predict the input flux, and so it is useless at high hght 
levels, as seen from the misfit function (Fig. (HJ. 

For the PP thresholding strategy, the ENF tends to- 
wards \/2 as the hght level increases, and above about 
twenty photons per pixel we see that there is little advan- 
tage in thresholding, with the raw L3CCD output giving 
similar noise performance. However, thresholding does lead 
to a significant improvement in noise performance, reducing 
the ENF up to twenty photons per pixel. 

We also see that the misfit function remains low when 
we apply these thresholding schemes (Fig. |HJ using a com- 
bination of the PC and PP thresholding strategies. The PC 
thresholding strategy only performs well for light levels up 
to about 0.5 photons per pixel, and should not be used above 
1 photon per pixel. 



3.5 Photometric correction 

Fig. 13 shows that a nonlinearity correction for flux will be 
required during post-processing for hght levels greater than 
about 0.1 photons per pixel when using a PC thresholding 
strategy, and at hght levels less than about 20 photons per 
pixel when using a PP thresholding strategy. We can cor- 
rect the flux for a PP thresholding strategy approximately 
according to 



(6) 



1 -f 0.7exp(-/oet/3) 



where last is the result of the initial thresholding and sum- 
mation process. The result of this correction is shown in 
Fig. 13 Similarly, the flux for a PC thresholding strategy can 
be corrected according to 

-^corrected = — lu (1 — /est) (7) 

though the error in corrected fiux becomes large as /est — > 1. 



3.6 Lucky imaging 

As an example of where thresholding techniq ues can be used, 
we c onsider the Lucky Imaging technique iBaldwin et alJ 
I2OOII) , as used in astronomy to overcome atmospheric effects 
on medium sized telescopes. Snapshot images are taken with 
very short exposure times (of order 10-30ms). Many such 
images are taken, and the images with least atmospheric 
smearing are kept and added together after centroiding. An 
L3CCD is required since the light levels will be very low. 

We can threshold each data frame using both the PC 
and PP thresholding strategies and summing with previ- 
ous frames immediately (in parallel since these thresholding 
schemes are applied after detection). At the end of the ob- 
serving run, these two images can then be combined after 
applying the nonlinearity correction for flux, depending on 
whether the signal on each pixel is low or high. 



3.7 Sources of errors 

Apart from when using the photon counting (single thresh- 
olding) strategy, we require knowledge of t he mean gain. 
This can be controlled to about 1 per cent iMackav et alJ 
12001!) which, with gains of order 1000, requires millivolt sta- 
bility for the clock-high L3CCD electrode voltage. Errors 
here will only have a small effect on our estimates at low 
light levels and at higher light levels, the fiux estimate error 
will be proportional to the error in gain. 



4 CONCLUSIONS 

We have studied the stochastic gain process of L3CCDs and 
characterised the noise added by this multiplication process. 
For different signal level regimes, we have investigated the 
best processing strategy for the L3CCD output with regard 
to estimating the true photon input. 
In summary, we find that: 

(i) Single thresholding of the L3CCD output can be ap- 
plied accurately for photon rates up to about 0.5 photons 
per pixel per readout with a small nonlinearity correction. 

(ii) Multiple thresholding (binning) of the L3CCD output 
can be applied at any light level, and is most advantageous 
with light levels between 0.5-20 photons per pixel per read- 
out. This can reduce the excess noise factor introduced by 
the multiplication process from \/2 to 1.1 at light levels of 
about one photon per pixel. 

(iii) Using a combination of single and multiple thresh- 
olding strategies leads to further improvement, decreasing 
the excess noise factor to unity at lower light levels. 

(iv) If the gain is not known accurately, threshold bound- 
aries will be chosen wrongly. However, the gain can be con- 
trolled to about 1 per cent, so this effect is small. 

Our recommendation when using L3CCDs at low light 
levels (less than 0.5 photons per pixel per readout) is that 
a single threshold processing strategy on the L3CCD out- 
put should be used. At higher light levels up to about 20 
photons per pixel, threshold boundaries placed with the PP 
thresholding strategy should be used. A nonlinearity correc- 
tion should then be used, leading to correct flux estimation 
and an improvement in SNR performance from •v/n/2 to 
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y/n/l.l in the best case, without requiring an initial esti- 
mate of the mean hght level. At high light levels greater 
than about 20 photons per pixel, using the raw output does 
not lead to worse SNR performance than that obtained with 
other thresholding techniques. 

Since LSCCDs provide the best SNR performance at 
low light levels using a single threshold, we recommend that 
if possible, they are always used in this regime, increasing 
the frame rate if necessary to keep the number of photons 
per pixel low (< 0.5). A new controller being developed for 
LSCCDs will allow pixel rates of up to SOMHz, allowing the 
potential of LSCCDs to be maximized. 
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Pn(c) ^ Pn-l{x) * Pl{x) = "^P„-l{x) X Pl{c ~ x) (A3) 

x — n 

where * represents convolution. This gives: 

{x -2){x+ 1) exp{-xg''^) 



Psix) = 



P4(x) = 



Pb{x) = 



332! 

E'^ — 7x — 6) exp{—xg^^) 



(A4) 



2x-^ - 13x^ - 278x + 936) exp{~xg'^) 



where Pn{x) is valid from x > n (and is in fact zero at 
X — n — 1). We can see that for moderately large x (and x 
usually will be large over most of the distribution, since the 
gain is large) we can simplify to give a general probability 
distribution: 



Pnix) = 



x" ^ exp{-xg ^) 



(A5) 



5" (n-1)! 

which has an expectation ng and variance ng^, and fits the 
actual distribution almost perfectly. Variations on this are 
possible, for example taking more care at small x, though 
since the overall differences are small, the simplified version 
is used here. This approximation is not valid for large n, and 
a Gaussian distribution with the same mean and variance 
should be used instead. 

This paper has been typeset from a Tp^X/ ffl^jX file prepared 
by the author. 



APPENDIX A: PROBABILITY DISTRIBUTION 

iMatsuo et alJ il985l) give the output probability distribution 
for a electron multiplication device (e.g. an LSCCD), with 
r elements in the multiplication register, and a probability 
P of producing an extra electron at each stage, for a single 
photon input, as 

Pr{x) = {l-P)pr-l{x) 

Pr{0) =0, r > 1 

po{x) = Si^x, r > 1 

If r is large and P is small, we find that this can be 
approximated to an exponential distribution: 

Pi{x) ^ g-^ exp{-xg-^) (Al) 

where -Pi (a;) signifies the probability of an output x for a 
single photon input, with mean gain g. This gives < x >= 
= g, as expected. 
To generate the probability distribution for two input 
electrons, we simply take the convolution of Eq. lAll with 
itself. This gives: 

c 

P2{c) = exp(-a:g~^) exp(-(e - a;)g~^) 

= g-^{x~l)exp{-xg-^) (A2) 

Likewise, the probability distributions for larger input 
electron counts can be derived using 



